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The low-temperature magnetic excitations of the two-dimensional spin-5/2 square-lattice Heisen- 
berg antiferromagnet Rb2MnF4 have been probed using pulsed inelastic neutron scattering. In 
addition to dominant sharp peaks identified with one-magnon excitations, a relatively weak contin- 
uum scattering is also observed at higher energies. This is attributed to neutron scattering by pairs 
of magnons and the observed intensities are consistent with predictions of spin wave theory. 



PACS numbers: 75.40.Gb, 75.30.Ds, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 



Rb2MnF4 is a near-ideal two-dimensional spin- 
5/2 Heisenberg antiferromagnet on a square lat- 
tice (2DHASL) and has been extensively studied 
experimentallji2i&4 to test theoretical predictions for 
this canonical 2D quantum Hamiltonian. Neutron scat- 
tering experiments observe sharp one-magnon excita- 
tions at low temperatures which could be well described 
by linear spin wave theory^ However, spin wave the- 
ory predicts that neutrons can also scatter from pairs 
of magnons leading to a broad scattering continuum at 
higher energies. Such processes are longitudinally polar- 
ized with respect to the ordered spin direction and cor- 
respond physically to the simultaneous creation of two 
magnons with opposite spin S z = ±1 (total AS Z = pro- 
cess). The overall scattering weight of the two-magnon 
continuum relative to one-magnon processes is related to 
the relative strength of the zero-point longitudinal quan- 
tum fluctuations in the ground state, which reduce the 
amount of ordered spin moment by AS compared to the 
full spin value S. Very sensitive experiments are required 
to search for two-magnon continuum scattering since for 
spin-5/2 the two-magnon scattering intensity integrated 
over energy and wavevector is expected to be only of the 
order 6% of the integrated one-magnon intensity. 

Another probe of two-magnon processes is light or Ra- 
man scattering and experiments on several 2DHAFSL 
systems have been made. Measurements in the S=l 
system K^NiF/f^ were found to be in excellent agree- 



ment with calculations based on interacting spin wave 
theory. Similar experiments on the cuprates (S=l/2) 
have, however, shown a clear discrepancy between the 
observed lineshape and calculations based on spin wave 
theory for pure 2DHAFSLii The discrepancy is believed 
to arise from the presence of the four-spin cyclic exchange 
terms which are also used to explain the observed spin 
wave dispersion along the antifcrromagnetic zone bound- 
ary in La2Cu04 ref. llCI 

Neutron scattering (unlike Raman scattering which has 
inherent momentum constraints) can in principle access 
the full wavevector and energy dependence of the two- 
magnon scattering. Raman and neutron scattering offer 
complementary information in that neutron intensities 
arc related to two-operator correlation functions, whereas 
Raman is related to four-operator terms. 

So far, two-magnon continuum scattering has been ob- 
served in the 3D material CoFzii*^, but little has been 
done to quantitatively measure the two-magnon scatter- 
ing using neutrons in 2DHAFSL systems, the only other 
experiments we are aware of are on the spin-1/2 material 
Cu(DCOO) 2 .4D 2 (CFTD)pi2i The large quantum cor- 
rections in 5= 1/2 systems make it important to measure 
the two-magnon scattering in Rb 2 MnF4 with 5 = 5/2 
where quantum fluctuations are expected to be smaller 
and to test the extent to which spin wave theory can 
describe the results. 

The purpose of this paper is threefold. Firstly we 
wished to measure the spin waves with the MAPS spec- 
trometer at ISIS to investigate whether high quality mea- 
surements of the dispersion relation could be made with 
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a neutron time-of-flight spectrometer as compared with 
earlier triple-axis measurements^ Secondly we wished to 
study the spin wave energy along the antifcrromagnctic 
zone boundary to search for evidence of four-spin inter- 
actions, second neighbour exchange interactions or quan- 
tum corrections to linear spin wave theory. Thirdly we 
wished to look for the two-magnon scattering continuum 
predicted by spin wave theory but lacking experimental 
evidence. 

The rest of the paper is organized as follows. In Sec. 
ITT1 we review the spin wave theory predictions for one- 
and two-magnon scattering. Details of the experiment 
are given in Sec. IIIII followed in Sec. IIVI bv results for 
the main dispersion relation with particular attention 
to wavevectors on the antiferromagnetic zone boundary. 
Sec.EJreports on the observation of two-magnon scatter- 
ing. The results arc summarized in a final Sec IVII 



II. DYNAMICAL CORRELATIONS 

Rb 2 MnF4 crystallizes in the tetragonal K 2 NiF 4 struc- 
ture with space group I4/mmm and lattice parameters a 
= b = 4.215A, c = 13.77A. The magnetic ions are Mn 2+ 
with a spin-only moment of S = 5/2, arranged in a square 
lattice in the basal plane, with antifcrromagnctic su- 
perexchange couplings between nearest neighbours medi- 
ated by intervening F~ ions. The inter-plane coupling is a 
factor of 10 -4 of the intra-plane coupling^. This almost 
perfect two-dimensionality arises because MnF 2 mag- 
netic layers have a large separation along the c-axis filled 
by two non-magnetic RbF sheets, and furthermore the 
magnetic couplings along c are frustrated because each 
Mn 2+ ion is equidistant to four antifcrromagnctically- 
coupled spins in the layers below and above, further 
weakening the effect of the inter-plane couplings. 

Antiferromagnetic order with moments along the c-axis 
occurs below Tn = 38.4 K. Earlier studies^ proposed 
that the observed ordering at finite temperature can be 
well accounted for by a small anisotropy S z = 0.0048(10) 
in the Hamiltonian, ultimately originating from dipolar 
interactions, i.e. 



H = JJ21S?S? + SfSf, + (1 + S Z )S*SZ 

(u>) 



(1) 



where J ~ 0.6544(14) meV [ref. Q is the exchange energy 
for nearest-neighbour spins on the square lattice and {ii') 
indicates that each interacting spin pair is counted once 
in the summation. (x,y,z) are along the crystallographic 
(a,6,c) axes. For this Hamiltonian the magnon dispersion 
relation in linear spin wave theory is^ 
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where jq = cosn(Q k + Qh) cos7r(-Q fc + Q h ) and 
(QinQk,Qi) are components of the crystal momentum 



Q given in rlu units of (27r/a, 2n/b, 2w/c). Often a mul- 
tiplicative factor, Z c = (1 + 0.157/25'), is included in the 
dispersion relation to account for corrections to lowest or- 
der spin wave theory. For S = 5/2, Z c = 1.0314, and we 
neglect this correction as it can be readily incorporated 
into the exchange constant. 

Neutron scattering measures the dynamical correlation 
functions given by 



S aa (Q, ui) ~ 
2*hNJ-oo if 



^Q-{ r i— r i 



(3) 



where TV is the total number of spins and the sum runs 
over all sites i and i' in the lattice. One-magnon events 
occur in the spin correlations transverse to the ordered 
spin direction z. In the non-interacting spin wave ap- 
proximation the transverse correlations at T = K are 
given by 1 ^ 



S**(Q,lo) 



syy{Q,u) 

\{S - AS)(u Q + v Q ) 2 8{hw - fiwg) (4) 



cosh( 



sinh^, tanh20 = — jj 



1Q 



where u Q - „ w - - (T+57)' 

Here AS" = S — (S z ) is the spin reduction due to zero- 
point fluctuations calculated as (1/N) J2q v q where the 
sum extends over Q's in the full Brillouin zone. AS = 
0.197 for the isotropic Heisenberg model and 0.167 for 
the anisotropy S z appropriate for Rb2MnF4. 

The finite spin reduction allows for the presence of lon- 
gitudinal fluctuations, which can be described in terms 
of two-magnon scattering events. In the non-interacting 
spin wave approximation the longitudinal correlations at 
T = K are±£ 

S ZZ (Q, ^elastic = ^ /(Ql,Q2)x 
Qi,Q 2 

&{Q - Qi + Q2) 5{hw - fnoQ 1 - fujjQ. 2 ) (5) 

where f(Qi,Q2) = \ {uq 1 vq 2 - UQ 2 VQ t ) 2 is the struc- 
ture factor for creating two magnons at wavevectors Q\ 
and Q 2 - In the summation, one of the two magnons 
(say Qi) is restricted to the first Brillouin zone. The 
above equation gives the inelastic part of the longitu- 
dinal correlations, the elastic part is simply the Bragg 
peak contribution (S — AS) 2 S(fku) S(Q — Qaf — t), 
where Qaf=(0.5,0.5) is the antiferromagnetic ordering 
wavevector and r is a vector of the reciprocal lattice. 

An understanding of how the scattering cross-section 
is distributed between the elastic, one and two-magnon 
channels can be obtained by comparing the integrated 
scattered intensities. The total intensity integrated over 
energy and averaged over the Brillouin zone is given by 
the well known sum rule J S(Q,u)dQd(huj) = S(S + 1). 
Similar expressions can be derived for the individual scat- 
tering components, by integrating over the expressions 
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for the elastic, one-magnon eq. J3J and two-magnon eq. 
components. The results are summarised in Table [I] 



component 


integrated intensity 


S(Q,u) 
S ZZ (Q, w) c i ast i c 

S" Z (Q, ^)inclastic 


S(S + l)= 8.75 
(S - AS) 2 = 5.443 
(S- AS)(2AS + 1) =3.112 
AS(AS + 1) = 0.195 



TABLE I: Total sum rules for the different components of the 
scattering, evaluated for S = 5/2 and AS = 0.167. 



The inelastic neutron scattering intensity from one and 
two-magnon excitations is proportional to 

(2- Pz )S xx (Q,u>) + Pz S zz (Q,Lu) inclastic (6) 

where p z = 1 — Q z is the polarization factor for lon- 
gitudinal scattering, Q z is the directional cosine of the 
wavevector Q with respect to the c-axis. The proportion- 
ality factor between eq. @) and neutron scattering inten- 
sity includes the magnetic form factor squared F(Q) 2 of 
Mn 2+ ions. 



III. EXPERIMENTAL DETAILS 

The magnetic scattering was measured from a 13. 4g 
single crystal of Rb 2 MnF4 using the MAPS spectrometer 
at the ISIS pulsed neutron source in the UK. The sample 
was enclosed in an aluminium can containing helium ex- 
change gas and measurements reported here were made 
at the base temperature of 9.5 K. MAPS is a direct ge- 
ometry time-of- flight instrument, equipped with a 16 m 2 
array of position sensitive detectors, divided into nearly 
4 x 10 4 separate detector elements. This allows collection 
of the inelastic scattering intensity in a highly pixelated 
3D volume in the 4D (Qh, Qk, Qi, oj) space, from which 
one can extract the intensity plot in a certain plane or 
along a certain direction as illustrated in Fig. and 
c). We describe the data in terms of the two in-plane 
wavevectors (Qh, Qk) and energy fku as the magnetism is 
two-dimensional and the inter-layer component Qi only 
enters the magnetic scattering through the magnetic form 
factor F(Q) and the polarization factor p z with respect 
to the ordered spin direction, see eq. iJBJ. 

An incident neutron energy of 24.92 meV was selected 
to map the inelastic scattering over the whole dynamic 
range of one- and two-magnon scattering processes which 
extended up to 13 meV. A Fermi chopper spinning at 300 
Hz gave an energy resolution of 0.75± O.OlmeV (FWHM) 
on the elastic line. Measurements were made with the 
two-dimensional magnetic layers arranged perpendicular 
to the incident neutron beam (orientation 1, c || fc^, a 
horizontal) and tilted by an angle tp = 45° with respect 
to the incident beam direction (orientation 2, (c, fej) = 0) 
to collect complementary data with different longitudinal 



polarization p z at the same two-dimensional wavevector 
and energy (Qh,Qk,u). Typical counting times for one 
crystal orientations were 20 hours at an average proton 
current of 170 /iA. To increase statistics the data was 
folded along symmetry-equivalent axes as illustrated in 
Fig. |21 four-fold in the c-axis along ki setup, and two- 
fold in the rotated configuration. The low angle detector 
bank 29 < 30° provided coverage over most of the first 
Brillouin zone. 



IV. SPIN WAVE DISPERSION RELATION 

The magnon dispersion relation was determined from 
fits to energy scans at constant wavevector (Qh,Qk) and 
typical data is shown in Fig. [21 The scattering is domi- 
nated by a sharp one-magnon excitation and solid lines 
show fits to eq. I@J), where the delta function S(Huj — Hujq) 
is replaced by a resolution broadened Gaussian peak. 
The extracted spin wave dispersion along symmetry di- 
rections in the Brillouin zone is shown in Fig. 01 Data at 
the lowest energies is limited because of the difficulty in 
resolving the one-magnon peak from the elastic incoher- 
ent scattering. The solid line shows a fit to the dispersion 
relation in eq. with a fitted exchange J = 0.648±0.003 
meV and fixed anisotropy 8 Z = 0.0048, in agreement with 
the previous estimates of J = 0.6544 ± 0.0014 meV ob- 
tained from triple-axis neutron measurements^. 

Linear spin wave theory, eq. @ , predicts no dispersion 
along the antiferromagnetic zone boundary (central panel 
in Fig. 0} therefore an observed dispersion along this di- 
rection can give information about failures of the lin- 
ear spin wave approximation or extra interactions in the 
Hamiltonian beyond the nearest-neighbour exchange. In 
the spin- 1/2 organic material CFTDi&, the zone bound- 
ary energy was observed to decrease from (0.25,0.25) to 
(0.5,0) by 6 ± 1%, and this agreed with computational 
work on the 5=1/2 2DHAFSL using the Ising-limit 
expansion that predicts a dispersion of 7%i£. The zone- 
boundary dispersion was therefore attributed to quantum 
corrections to linear spin wave theory. A spin wave calcu- 
lation extended to order l/S 2 predicts the same sign but 
a much smaller magnitude (2%) of this dispersion^ than 
the series results, suggesting that higher order corrections 
would need to be considered to obtain a fully satisfactory 
theory. In La2Cu04 the observed zone-boundary disper- 
sion had the opposite sign, increasing from (0.25,0.25) to 
(0.5,0) by 13% i£. This effect was attributed to higher- 
order spin exchange terms in the Hamiltonian. La2CuC>4 
is a Mott insulator and a more appropriate description of 
the electronic states is in terms of a Hubbard model at 
half-filling, characterized by a kinetic energy gain from 
hopping, t, and potential energy cost, U, for two elec- 
trons occupying the same site. At small t electrons are 
localized, and perturbative expansion in t/U gives an 
effective Hamiltonian for the spin degrees of freedom. 
The first term in this expansion is the nearest-neighbour 
Heisenberg exchange, and the dominant next order term 
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FIG. 1: (Color online) (a) One-magnon dispersion surface as a function of two-dimensional wavevector (Qh,Qk) and energy 
ui (color shading is intensity in neutron scattering). Dashed lines in the basal plane and at maximum energy hui = 4JS mark 
the antiferromagnetic zone boundaries. The basal plane also shows constant-energy contours (solid lines), (b) Constant energy 
maps of the magnetic scattering at huj = 1, 3.5 and 6 meV obtained by taking slices from the 3D (Qh,Qk,u) neutron data, (c) 
Energy scan at constant wavevector (—0.5, 0) along the direction shown by the vertical rectangular column in (b) (cross-section 
of column AQh x AQ^ indicates the wavevector region over which intensity points were averaged). 




orientation 1 slice 



orientation 1 slice after folding 



a o 



a ° 




FIG. 2: (Color online) Folding of data along symmetry- 
equivalent directions: (a) four-fold symmetry in the c || ki 
setup (b) two-fold symmetry in the rotated configuration. 
Data corresponds to intensity at an energy 5.5 ± 0.5 meV. 



is a cyclic exchange coupling 4 spins at the corners of 
each square plaqucttc. Such a ring exchange term corre- 
sponding to U/t ~ 6 was used to describe the high-energy 



magnon dispersion^ Raman scattering^ and infrared ab- 
sorption experiments^ in La2Cu04. 

The spin wave energies along the antiferromagnetic 
zone-boundary contour in Rb2MnF4 are shown in a 
blown-out scale in Fig. [5] The dispersion along the 
Qh + Qk = 0.5 direction is very small, of order 1 ± 0.5% 
between (0.5,0) and (0.25,0.25). In extracting peak po- 
sitions we carefully considered the effects of wavevector 
averaging over a box of finite size AQh x AQk around the 
nominal (Qh, Qk) values. Since the dispersion surface has 
a maximum on the zone boundary, the effect of a finite 
wavevector averaging is to produce an apparently lower 
peak energy. The effect is more pronounced around the 
(0.25,0.25) point since there the dispersion surface has 
only a one-dimensional maximum, whereas the corner 
point (0.5,0) is a local maximum along two directions 
in the plane, see Fig. ^t). This effect is illustrated in 
Fig.HJlby taking cuts of different sizes AQh x AQk over a 
simulated data set for the one-magnon scattering cross- 
section, eq. (0J. A wavevector box of size 0.04 x 0.04 was 
chosen as a balance between a minimal apparent peak 
shift (< 0.1%) and sufficient data pixels in the box for 
sufficient statistics to allow quantitative analysis and the 
final results are shown in Fig. [SJ Data collected under 
two sample orientations show a small dispersion with the 
higher energy at (0.25,0.25). The magnitude of the dis- 
persion is small, 1 ± 0.5% of the zone boundary energy, 
a value close to the limit of the experimental accuracy, 
which may explain why the two data sets are not exactly 
overlapping. 
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FIG. 3: Energy scans at constant (Qh,Qk) wavevector. Solid 
lines are fits to Gaussian peaks. 
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FIG. 4: Dispersion relation along symmetry directions in the 
Brillouin zone (bold lines in Fig. 0b)). Solid line shows a fit 
to eq. © with J = 0.648 meV. 



We note that quantum corrections to linear spin wave 
theory to order 1/S 2 would predict a dispersion of the 
same sign but nearly an order of magnitude smaller, 0.005 
meV— . Earlier triple-axis measurements also observed a 
dispersion similar to the one in the present experiments 
and Cowley et alA proposed that the origin was a next- 
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FIG. 5: Magnon energy for wavevectors along the antiferro- 
magnetic zone boundary obtained using narrow (0.04 x 0.04) 
cuts from data collected under two different experimental se- 
tups. Solid line is a fit to the spin wave dispersion relation 
including next-nearest neighbour coupling, eq. (0). 
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FIG. 6: Energy of one-magnon peak extracted from cuts 
through a simulated data set for eqs. 12141 . Using a large 
wavevector averaging range AQh X AQk = 0.1 X 0.1 around 
the nominal (Qh,Qk) position introduces an apparent dis- 
persion of 0.9% (filled circles), whereas the smaller cut size 
0.04 x 0.04 (open symbols) reduces this to less than 0.1%. 



nearest neighbour antifcrromagnetic exchange J' along 
the square diagonals. The dispersion relation in this case 
becomes 



fkt), 



Q 



US 



2J 



7q 



1/2 



(7) 



where 7q = cos7r(Q fc + Q h ) cosn(-Q k + Qh), 7q = 
cos[27r((3? l + Qfe)]+cos[27r(-Q/ l + (3 A; )]-2. We have fitted 
this expression to the dispersion along the zone bound- 
ary in Fig.0and the results are J = 0.657 ± 0.002 meV, 
J' = 0.006 ± 0.003 meV, compared to J = 0.673 ± 0.028 
meV, J' = 0.012 ±0.002 meV obtained previously! The 
difference in J arises because the frequencies at the zone 
boundary are lower in our measurements due possibly to 
small errors in the absolute energy calibration in one of 
the two experiments. 



6 



In short the present experiments observe a defi- 
nite change in energy along the antiferromagnetic zone 
boundary. The effect is small and the data barely pro- 
duces a reliable estimate of the effect. Nevertheless, our 
results do suggest that its most probable origin is from 
the next nearest neighbour exchange constants. We shall 
however neglect this small effect for most of the remain- 
der of this paper. 



V. TWO-MAGNON CONTINUUM 
SCATTERING 

To estimate the two-magnon scattering intensity pre- 
dicted by spin wave theory we evaluated eq. JHJ) numeri- 
cally by summation over a grid of finely spaced Q\ points 
over one Brillouin zone, setting Q2 = Q — Qi and replac- 
ing the delta function in energy with an area-normalized 
narrow Gaussian. Fig. UJi) shows an overview plot of the 
calculated two-magnon cross-section as a function of en- 
ergy and wavevector along symmetry directions the Bril- 
louin zone. Two-magnon scattering occurs in the form 
of a continuum at higher energies above the one-magnon 
dispersion relation. The small energy separation E sap (Q) 
between one- and two-magnon excitations is a conse- 
quence of the small uniaxial anisotropy S z in the Hamilto- 
nian eq. which opens a gap in the one-magnon spec- 
trum at the zone centre ftioa = AJ S\/8 Z (2 + S z ), with 
two-magnon scattering starting at the higher energy of 
2 x fiwo- Generally the two-magnon intensity is strongest 
for low energies and wavevectors near the antiferromag- 
netic zone centre, but here is also where the one-magnon 
structure factor is largest. The clearest way to separate 
a scattering signal from one- and two-magnon processes 
is at energies above the one-magnon zone boundary, and 
Fig. EJd) shows a plot of the calculated two-magnon in- 
tensity distribution in the Brillouin zone at an energy 
8.75 ± 1.25 meV, above the one-magnon cut-off. 

Calculated lineshapes at representative points in the 
Brillouin zone are shown in Fig. [51 One noticeable feature 
is the appearance of strong singular peaks, which become 
more prominent upon increasing the numerical accuracy 
in evaluating eq. JSJl. The singularities are a result of 
divergencies in the two-magnon density of states obtained 
by putting / = 1 in eq. (|SJ), i.e. 
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Plots of D(Q,u>) are shown by dotted lines in Fig. |H1 
The singularities at the high energy boundary present in 
D{Q,uS) do not show up in S ZZ (Q, cj)i ne i a stic- This is be- 
cause the structure factor of those processes in neutron 
scattering cancels as both magnons are on the antiferro- 
magnetic zone boundary contour where VQ t = vq 2 = 0. 
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FIG. 7: (Color online) (a) Two-magnon scattering intensity 
S ZZ (Q, w)i ne iastic as a function of energy and wavevector along 
symmetry directions in the Brillouin zone (thick solid curve 
in (b)). Colour represents intensity. The solid line shows the 
one-magnon dispersion relation Hluq from eq.(|5J. (b) Two- 
magnon intensity as a function of wavevector at constant en- 
ergy huj — 8.75 ± 1.25 meV. Lineshapes at positions labelled 
A-F are shown in Fig. |H1 



It is interesting to consider whether any of the sin- 
gularities are a true feature of the two-magnon neutron 
scattering or whether they are a consequence of using 
non-interacting spin wave theory. Canali and WallinSS 
have included first order spin wave interactions in calcu- 
lating S zz (Q,ui) and their results show that the singu- 
larity peaks remain. However, any treatment of interac- 
tions within spin wave theory is perturbative and so it is 
possible higher order interactions may be still be impor- 
tant when considering the singularities. Experimentally, 
resolution effects would make the singularities very diffi- 
cult to observe as shown in Fig. [S] where the two-magnon 
lineshapes are convoluted with the resolution in our ex- 
periment and the results (dashed lines) show that any 
singularity would probably not be visible in our mea- 
surements. 

To test for the presence of two-magnon scattering in 
the data, we extracted energy cuts at fixed wavevectors 
(Qh,Qk)- Because of the inherently weak intensity of 
the two-magnon cross-section we chose a wavevector av- 
eraging range 0.1 x 0.1 to have enough data pixels for 
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FIG. 8: Calculated two-magnon scattering lineshapes at fixed 
wavevectors, indicated by labels A-F in Fig. |7| Solid line 
is S ZZ (Q, w)ineiastic and dashed lines show the effects of the 
instrumental resolution. Dotted lines show the two-magnon 
density of states, D(Q,lj) eq. JHJl, divided by a factor of 10. 



quantitative analysis. A typical scan near the antiferro- 
magnetic zone boundary at (0.5,0.15) is shown in Fig. 

The inelastic scattering is dominated by a sharp, 
one-magnon peak centred at 6.40 ± 0.02 meV, and ad- 
ditional much weaker scattering is observed in the form 
of a high energy continuum tail extending to at least 9 
meV [see Fig. EJi)], much higher that the one-magnon 
zone boundary energy. 

Our approach is to fit the data to a lineshape contain- 
ing both one and two-magnon components, with their 
relative intensity fixed by theory, eqs. 10} , J5j and JfjJ). 
The effects of resolution broadening are also included as 
discussed below. The effects of the finite cut size were 
included by averaging the predicted intensity over the 
finite wavevector size AQh x AQfc of the cut. The re- 
sulting profile was then convolved with the energy res- 
olution of the spectrometer. This was determined from 
the observed lineshape of the quasi-elastic peak in Fig. 
|5Jd). This showed a slightly asymmetric tail at lower ener- 
gies (due to the asymmetric neutron pulse shape) and the 
whole profile could be well parameterized by a sum of two 
Gaussian peaks, one off-centred on the low-energy side, 



see Fig. 0;). Such a weakly-asymmetric lineshape also 
provided a good description of the observed one-magnon 
peak lineshape as shown in Fig. |3Ji) . The relative posi- 
tions, intensities and widths of those two Gaussians are 
fixed while the fitted parameters are the overall width of 
the one-magnon peak (to parameterize the variation of 
the energy resolution with energy transfer) and an over- 
all scale factor. Fig. shows that the whole observed 
scattering lineshape including the high-energy tail can 
be well described when the two-magnon cross-section is 
included. 

a) Q=(0.5,0.15) - full cut lineshape 




Energy (meV) 

d) one-magnon tails - 
fit without two-magnon 
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FIG. 9: Stages of the fitting procedure, a) Energy scan at 
(0.5,0.15). b)-c) Low-energy asymmetric lineshape of the in- 
coherent scattering is well parameterized by a two Gaussian 
lineshape. d) Fitting the same resolution profile to the one- 
magnon peak centred at 6.40 meV, the extra scattering on the 
high-energy tail can not be accounted for by resolution effects, 
e) Fits to the combined one- and two-magnon scattering eq. 
© including resolution effects. 

Fig. 1101 shows the overall comparison between the ob- 
served intensity in the high-energy tails and that ex- 
pected from two-magnon scattering. The plot corre- 
sponds to intensity averaged over the energy range 7.5 
to 10 meV, much higher than that of the one-magnon 
zone boundary, where only two-magnon processes are ex- 
pected to contribute. There is good overall agreement 



8 



between data and calculations, which include the longi- 
tudinal polarization factor p z expected for two-magnon 
processes eq ©■ 




Q 2 4 6 8 2 4 0.6 8 




0.2 4 0.6 0.8 0.2 0.4 0.6 0.8 



Q h Q h 

FIG. 10: (Color online) Comparison between data and pre- 
dicted two-magnon scattering at an energy hui = 8.75 ± 1.25 
meV, much higher that the one-magnon zone boundary en- 
ergy, i) and iii) show data in two different experimental config- 
urations and ii)and iv) show the simulated two-magnon scat- 
tering intensity including the polarization factor, magnetic 
form factor, an overall scale factor and a flat background to 
compare with data. Boxes labelled a)-f) show location of en- 
ergy scans plotted in Fig. 1111 The dashed square box is the 
antiferromagnetic zone boundary. 

A number of representative energy scans extracted 
from the data near the antiferromagnetic zone bound- 
ary where the two-magnon contribution can most easily 
be singled out are shown in Fig. ^2 Solid lines show the 
results of the fitting procedure described above to the 
combined one-and two-magnon scattering lincshapes eq. 
(jSJ) and good agreement is observed throughout. In those 
fits the overall scale factor was allowed to vary for each 
scan, however very similar fits are obtained if this scale 
factor is fixed to a common average value for all scans 
and for clarity we only show the two-magnon component 
determined in this way (dashed lines), practically indis- 
tinguishable from the results of the free fits. As an inde- 
pendent consistency check we have converted the neutron 
intensities into absolute units of barn meV -1 sr" 1 per 
spin normalizing by the sample mass and by intensities 
measured with a vanadium standard; the calculated two- 
magnon neutron scattering intensities in absolute units 
are very close (within 10%) to the results of the fits, this 
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agreement giving further support to the identification of 
the continuum intensities with two-magnon scattering. 

The scans shown in Fig. ^2 include points where the 
intensity from two-magnon scattering is predicted to be 
low, as well points where it is predicted to be high, this 
modulation of intensity being mainly due to the polar- 
ization factor p z in eq. 10 • For example scans in Fig. 
II lb Vf ) have rather weak high-energy tails and cannot 
be taken in isolation to provide evidence of two-magnon 
scattering, but are significant when taken in the context 
of all the scans shown. The scans in Fig. Illb'Vb') show 
the strongest high-energy signal, since they correspond 
to the largest longitudinal polarization factor appropriate 
for two-magnon scattering (and lowest polarization factor 
for transverse one-magnon scattering). Another illustra- 
tion of the polarization effect is provided by comparing 
the data in Figs.lTTb') and d') collected at equivalent po- 
sitions in the Brillouin zone: using the same two-magnon 
intensity in both scans without adjusting for the change 
in polarization factor results in a large disagreement (dot- 
ted line in Fig. Illb ') with the observed continuum scat- 
tering intensity. From this analysis we conclude that the 
observed continuum scattering at high energies is consis- 
tent both in magnitude and polarization with scattering 
expected from two-magnon processes, neglecting any in- 
teractions between the magnons. 



VI. CONCLUSIONS 

We have conducted a detailed investigation of the low 
temperature dynamical properties of the square-lattice 
spin-5/2 Hciscnberg antifcrromagnct Rb2MnF4. The 
spin wave dispersion was measured and a small varia- 
tion in energy along the antiferromagnetic zone boundary 
with found. The energy change along the zone bound- 
ary was too large to be explained by quantum correc- 
tions to linear spin wave theory, and could be the re- 
sult of weak next-nearest neighbour exchange interac- 
tions (J' I J = 1 ± 0.5%). Furthermore, a low intensity 
signal was observed around the high energy tail of the 
one-magnon peaks. The lineshapc and intensity variation 
of this signal provides good evidence that it is the result 
of scattering by pairs of non-interacting spin waves (two- 
magnon scattering). We conclude that although spin 
wave interactions are important in describing the shape 
of the two-magnon Raman scattering, they are much less 
important for two-magnon neutron scattering. 

This research was supported by UK EPSRC (TH, RC, 
RAC, DAT) and by US NSF through grant number 
DMR-0134377 (RLL). We are also grateful to the staff 
at ISIS for their help and support and to Bella Lake and 
Robert Birgeneau for useful discussions. 
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FIG. 11: (Color online) Energy scans analysed with a combined one and two-magnon cross-section eq. © (solid lines). Dashed 
lines shows the two-magnon component, with a common overall scale factor for all scans, as used in the simulation in Fig. 
1101 Scan labels a-f ) (with / for data in orientation 2) refer to locations in the Brillouin zone indicated in Fig. 1101 p z is the 
longitudinal polarization factor in the middle of the two-magnon continuum region at 8.75 meV. Dotted line in b') illustrates the 
type of disagreement obtained if one uses the two-magnon intensity observed at the equivalent position d') without adjusting 
for the change in polarization factor. 
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